arXiv:astro-ph/0503339vl 16 Mar 2005 


Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 


Spectroscopy of the type la supernova SN 2002er: days —11 

to +215 

R. Kotak 1 , W.RS. Meikle 1 , G. Pignata 2 , M. Stehle 3 , S.J. Smartt 4 , S. Benetti 5 , W. Hillebrandt 3 , D.J. Lennon', 

P.A. Mazzali 3,6 , F. Patat 2 , and M. Turatto 5 

1 Astrophysics Group, Imperial College London, Blackett Laboratory, Prince Consort Road, London, SW7 2AZ, U.K. 

2 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei Miinchen, Germany 

3 Max-Planck-Institut fur Astrophysik, P.O. Box 1317 D-85741 Garching, Germany 

4 Department of Physics and Astronomy, Queen’s University Belfast, Belfast, BT7 INN, Northern Ireland, U.K. 

5 Osservatorio Astronomico di Padova, vicolo dell’Osservatorio 5,1-35122 Padova, Italy 

6 Osservatorio Astronomico di Trieste, Via Tiepolo, 11,1-34131 Trieste, Italy 

' Isaac Newton Group of Telescopes, Apartado de Correos 321, 38700 Santa Cruz de la Palma, Canary Islands, Spain 
Received; Accepted; 

Abstract. We present an extensive set of optical spectroscopy of the nearby type la supernova, SN 2002er, with 24 epochs 
spanning — 11 to +34 days. Its spectral evolution is fairly typical of a type la supernova although it suffers high extinction. 
Nevertheless, there are differences in the spectral evolution when compared to coeval spectra of other normal type la supernova 
with comparable decline-rate parameters. Modelling of the photospheric phase spectra using a homogeneous abundance distri¬ 
bution in the atmosphere provides a fair match to the observations, but only by pushing the adopted distance and risetime close 
to the observational limits. Future improvements here will require models with a more realtistic stratified abundance distribu¬ 
tion. From simple modelling of a nebular spectrum obtained at +215 d, we infer a 56 Ni mass of 0.69Mq, consistent with that 
derived from the light curve. 

Key words, stars: supernovae: optical spectroscopy stars: supernovae: individual: SN 2002er 


1. Introduction 

Type la supernovae (SNe la) are currently the most accurate 
high-z cosmological distance indicator. Through the light curve 
shape corrections, accuracies of better than 10% in distance are 
claimed (e.g. Phillips et al. 1999; Riess et al. 1998). However, 
a basic assumption is that the correlations hold at all universal 
epochs. Yet, even for the best nearby SN la sample it appears 
that the various incarnations of the light curve correction meth¬ 
ods may not be wholly consistent with each other (Drell et al. 
2000; Leibundgut 2000). 

Although there is a general consensus that type la su¬ 
pernovae result from the thermonuclear fusion of C-O white 
dwarfs, observational support for the numerous progenitor 
channels that have been proposed (e.g. Branch et al. 1995) is 
only just beginning to be gathered (e.g. Ruiz-Lapuente et al. 
2004; Hamuy et al. 2003; Kotak et al. 2004). In contrast, there 
is an increasing body of evidence for the observational diver¬ 
sity (e.g. Benetti et al. 2004b) among type la supernovae per¬ 
haps suggestive of diversity in the progenitor channels and/or 
explosion mechanism. 

In the favoured scenario leading to a type la explosion, 
the white dwarf is one component of a close binary, and ac- 
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cretes hydrogen from its companion star. Other scenarios in¬ 
clude the merging of binary white dwarfs, or the accretion of 
He rather than H (e.g. Hillebrandt & Niemeyer 2000). When 
the white dwarf reaches the Chandrasekhar limit, explosive 
carbon burning sets in. Burning to nuclear statistical equilib¬ 
rium ensues, yielding mostly radioactive 56 Ni. At lower den¬ 
sities, intermediate-mass nuclei, e.g. 28 Si, are produced. Thus 
the typical observed early-time spectra of SNe la are domi¬ 
nated by lines of both iron-group and intermediate mass el¬ 
ements. However, the exact nature of the explosion mecha¬ 
nism is uncertain. Competing models can be tested using de¬ 
tailed spectral observations since they are highly sensitive to 
the nature of the explosion. For example, much higher veloc¬ 
ities in intermediate-mass elements are observed in delayed- 
detonation model spectra than in those of deflagration mod¬ 
els (e.g. Khokhlov 1991). Although the optical region is spec¬ 
troscopically crowded and line-identification is difficult due to 
the additional complication of blending as a result of Doppler 
broadening (several thousand kms '), it is also a region rich in 
diagnostic potential. For instance, based on measurements and 
modelling of lines such as Sill A6355 A, CII A6578, A7231 A, 
and the Call H and K lines, Benetti et al. (2004a) argue that the 
combination of high-velocity intermediate mass elements, an 
unusually low temperature, and low abundance of carbon indi- 
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Fig. 1 . B band image of SN 2002er taken around maximum light. The 
dashes indicate the position of the supernova. 


Table 1 . Properties of SN 2002er and its host galaxy 

Host galaxy UGC 10743 

Host galaxy type Sa 1 

RA (2000) 17'*ll m 29 s .88 

Dec (2000) +07°59'44".8 

Recession velocity (km s~ 1 ) 2569 

Distance modulus* 32.9 ± 0.2 

(Ho = 71kms _1 Mpc _1 ) 

E (B - V)* 0.36 

Offset from nucleus 12".3W 4".7N 

Date of B maximum (MJD)* 52524.2 (Sep. 07, 2002) 

Magnitude at maximum* U=14.72, B=14.89, V=14.59, 

R=14.43,1=14.49 

Amis(B)* _ 1.33±0.04 _ 

Notes. * Taken from Pignata et al. (2004); f NED. 

days later by Smartt et al. (2004); comparison with SN 1994D 
suggested that it was discovered at an early phase (about —10 d 
before maximum light) and that it was significantly redder. A 
list of basic parameters is given in Table 1. 


cate that burning to Si penetrated to much higher layers than 
expected for ‘normal’ type la supernovae. This was interpreted 
as evidence for a delayed-detonation explosion. An alternative 
scenario could be additional nuclear burning in the so-called 
distibuted regime at low densities which would have a similar 
effect (Ropke & Hillebrandt 2004). 

To date, good photometric and spectroscopic coverage ex¬ 
tending from pre-maximum light to late times has been ob¬ 
tained for only a handful of supernovae. In order to better 
understand and quantify the various peculiarities that have 
been observed, high quality photometric and spectroscopic data 
from pre-maximum through to nebular epochs is a prereq¬ 
uisite, with detailed comparisons with the predictions of ex¬ 
plosion models being made at every step. To this end, the 
European Supernova Collaboration (ESC) 1 aims to obtain un¬ 
precedented photometric and spectroscopic coverage for ~ 10 
nearby (v rec < 4000 kms -1 ) type la supernovae. The poten¬ 
tial of high quality early-time spectroscopy is illustrated by 
the method of “abundance tomography” (Stehle et al. 2004) 
whereby the detailed abundance stratification is derived by 
computing synthetic spectra for over ten epochs, using spec¬ 
troscopy from the earliest epochs to constrain later ones. 

We have already obtained excellent data sets for about eight 
nearby type la SNe. This is the third in a series of observational 
papers published on behalf of the ESC. We present early-time 
optical spectroscopy spanning 24 epochs of the second ESC 
target, SN 2002er. In a companion paper, Pignata et al. (2004) 
discuss the photometric evolution of SN 2002er. 

1.1. SN 2002er 

SN 2002er in UGC 10743 was discovered on 2002 Aug. 23.2 
U.T. in NEAT images by Wood-Vasey et al. (2002). and was 
spectroscopically confirmed to be a type la supernova three 

1 http://www.mpa-garching.mpg.de/Ttn/ 


2. Data acquisition and reduction 

The photospheric-phase observations span —11 to +34 days. 
Somewhat unusually, the spectroscopic monitoring began sev¬ 
eral days before the photometric campaign. Optical spec¬ 
tra of SN 2002er were obtained from several different tele¬ 
scopes (Table 2). The observations were carried out in service 
or target-of-opportunity mode usually using the configuration 
available at the time. This meant that the set-up varied with 
epoch and even with observations performed using the same 
instrument (see Table 2). 

The reduction of the data was carried out using standard 
NOAO IRAF and FIGARO routines (Shortridge 2002). The 
CCD frames were debiased and corrected for distortion and 
cosmic rays. Where a sufficient number of Hatfield frames were 
available, the data were also flatfielded. Sky subtraction was 
carried out by fitting low-order polynomials to the background 
on either side of the supernova spectrum and optimally extract¬ 
ing one-dimensional spectra from the two-dimensional frames. 
The wavelength solution was derived using arc frames and was 
checked against the bright night-sky emission lines. Flux cali¬ 
bration was carried out with respect to spectrophotmetric stan¬ 
dards usually taken using an identical instrument configuration 
and reduced as described above. No attempt was made to cor¬ 
rect for the telluric absorption features between 6000 - 8000 A. 

To obtain full wavelength coverage in the optical region, 
an observation at any given epoch typically consisted of two 
or three spectra taken either with a single grating/grism at dif¬ 
ferent central wavelengths or using different gratings/grisms. 
Severe distortion between the blue and red parts of the optical 
spectra was apparent on several nights and variable conditions 
exacerbated the situation. To correct for this distortion we mul¬ 
tiplied the fluxed supernova spectra by the BVRI Bessell filter 
functions and integrated across the bands. A similar procedure 
was applied to the spectrum of Vega, and hence spectroscopic 
magnitudes for the supernova were derived. Comparison with 
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Table 3 . E(B — V) values for the SNe shown in Fig. 3 



3 

to 

N 

Reference 

SN 1994D 

0.06±0.02 

Patat et al. (1996) 

SN 1996X 

0.1±0.03 

Salvo et al. (2001) 

SN 1998 aq 


Branch et al. (2003) 

SN 1998bu 

0.32±0.04 

Hernandez et al. (2000) 

SN 1999 ee 

0.32 

Hamuy et al. (2002) 

SN 2002bo 

0.43±0.1 

Benetti et al. (2004a) 


Notes. The absence of narrow Na D absorption lines in SN 1998 aq 
(Ayani & Yamaoka 1998) implies little reddening. 

near-contemporaneous photometry (from Pignata et al. 2004) 
was then used to derive scaling factors for each band. A low- 
order polynomial fit to the scaling factors was obtained which 
was then applied to full supernova spectrum. We estimate the 
flux calibration to be accurate to about 10-15%. For the earliest 
epochs (—11 to —8 d) no photometry was available for com¬ 
parison. Since it is the behaviour at the earliest epochs that we 
seek to understand, we did not scale these spectra in any pre¬ 
conceived manner. 

3. Optical spectra at early times 

This section is subdivided as follows: we detail the early spec¬ 
tral evolution of SN 2002er in Sect. 3.1 followed by a discus¬ 
sion on the presence of carbon in Sect. 3.2. We then describe 
the evolution of the expansion velocities in Sect. 3.3. 

3.1. Spectral Evolution 

The spectral evolution of SN 2002er from day —11 to +34 is 
shown in Fig. 2. The rapid evolution at early times can be fol¬ 
lowed in detail given the excellent temporal coverage. The out¬ 
ermost material in the ejecta, composed mainly of intermedi¬ 
ate mass elements is responsible for producing the spectrum 
up to about maximum light. The spectra shown in Fig. 2 are 
typical of normal type la supernovae, characterised at early 
epochs (up to about day +6) by the deep P Cygni absorption 
near 6150 A due to Sill and the distinctive W-shaped SII fea¬ 
ture near 5400 A. Other prominent features at these phases are 
attributed to the Call H and K blend, Mgll A 4481 A and a mul¬ 
titude of blends due to Fell, Sill, and SII in the 4500 to 5000 A 
range. The Nal D lines are clearly visible suggesting that sub¬ 
stantial reddening due to intervening material. 

A comparison with other type la supernovae at four repre¬ 
sentative epochs is shown in Fig. 3. The decline-rate parame¬ 
ters for these supernovae are fairly typical with the exception of 
the slow-decliner SN 1999 ee which shows markedly different 
behaviour at early epochs, but has been included for compari¬ 
son. Although differences between the spectra are present at all 
epochs, these are most pronounced in the pre-maximum spec¬ 
tra (Fig. 3a). The obvious contrasts are in the shape and depth 
of the prominent features and in the amount of structure in the 
spectrum. The latter is probably due in part to the differences in 
the velocities of the ejecta; this is especially true of SN 2002bo 
for which the higher velocities have probably smeared out the 


weaker features. However, the similarity between SN 2002er 
and SN 1994D, down to the weakest of features is striking. 
The strength of the temperature- sensitive Silll A4553.4568 A 
line suggests that SN 2002er and SN 1994D spectra are inter¬ 
mediate in temperature compared with the hotter SN 1998bu 
and SN 1998aq spectra and the cooler SN 2002bo spectrum. 
The W-shaped SII feature at A5454.5460 A is already apparent 
at —lid in the SN 2002er spectrum whereas it only became 
apparent at —2d in the spectrum of SN 1999ee (Hamuy et al. 
2002). 

By maximum light (Fig. 3b), the Call lines near 8500 A 
have become significantly stronger and the W-shaped SII 
A5454,5640A feature is very well-developed. The Call H & 
K feature (blended with the Sill A3858 A) is more similar to 
that of SN 2002bo than to the double-dipped feature seen in 
SNe 1994D, 1996X and 1998bu. Lentz et al. (2000) have sug¬ 
gested that the presence or absence of a split in this feature 
is related to the relative strength of the Sill and Ca features. 
The strength and width of the Sill A6355 A feature is compa¬ 
rable to that seen in SN 1994D and SN 1996X. The A4950A 
feature is stronger in SN 2002er at pre-maximum epochs than 
in the other supernovae shown in Fig. 3a; by maximum light, 
it has developed into a pronounced W-shaped feature. Unlike 
the W-shaped SII feature at ^5400A which is ubiquitous at 
maximum in all normal type la supernovae, only in SN 1996X 
does the A4950 A feature appear with comparable prominence 
to that of SN 2002er. The higher temperature of SN 2002er, 
as suggested by the Silll lines, lends credence to the sugges¬ 
tion that the A4950A feature is due to Felll (Mazzali et al. 
1995). Note that the presence of lines of high excitation is to 
be expected given the somewhat higher than normal 56 Ni mass 
inferred from the light curve by (0.7M. Pignata et al. 2004). 

Nugent et al. (1995) found a tight correlation between the 
ratio, 7?.(SiII) of the depths of the Sill A5972A and A6355A 
absorption troughs near maximum light and the luminosity of 
the supernova: the subluminous, faster decliners have larger 
values of 7?.(SiII). However, this correlation has been shown 
to break down for Atois(.B) < 1.2 (see Fig. 8 in Benetti et al. 
2004a). For SN 2002er we measure 0.23 for 77(SiIT) which is 
consistent with the Nugent/Benetti relation to within the errors. 

By about +20 d post maximum light, the Sill A6355 A fea¬ 
ture is still clearly visible (Fig. 3d). All the other supernovae, 
shown at comparable epochs, show kinks in the red wing of the 
Sill A6355 A feature due to Fell. SN 2002er, on the other hand, 
shows a rather smooth red wing. This suggests that the Sill 
feature is not confined to a narrow range in velocity and cor¬ 
roborates the finding by Pignata et al. (2004) from light curve 
modelling, that the ejecta are highly mixed. Notice also that the 
strength of the Call H&K feature is more similar to that of the 
slowest dec liner shown in Fig. 3, SN 1999ee. By about a month 
post-maximum light (Fig. 3e), the Sill feature is barely per¬ 
ceptible as features due to Fe dominate this wavelength range. 
Figs. 3d,e together suggest that the intermediate-element phase 
was of slightly longer duration in SN 2002er (i.e. slower transi¬ 
tion to Fe-peak elements) than for the other supernovae shown 
in Fig. 3. 

Interestingly, although the similarity of SN 2002er with 
SN 1994D and SN 1996X remains impressive throughout its 
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Table 2. Log of spectroscopic observations of SN 2002er 


Date 

M.J.D. 

Epoch* 

(d) 

Approx. Range 

(A) 

Telescope+Instrument+ 

Grating/Grism 

Flux 

Standard 

20020826 

52512.93 

-11.3 

3500-7200 

INT+IDS+R400V 

HD 19445 

20020828 

52514.87 

-9.3 

3600-7800 

INT+IDS+R632V 

Kopff 27 

20020829 

52515.86 

-8.3 

3300-8900 

INT+IDS+R632V 


20020830 

52516.84 

-7.4 

3600-9000 

A1.82+AFOSC+gm2,4 

BD+28 4211 

20020831 

52517.89 

-6.3 

3400-9900 

INT+IDS+R632V 

Kopff 27 

20020901 

52518.85 

-5.3 

3200-9100 

CA+CAFOS+B,R200 t 

BD+28 4211 

20020902 

52519.88 

-4.3 

3250-9300 

INT+IDS+R632V f 

Kopff 27 

20020903 

52520.85 

-3.3 

3200-10000 

CA+CAFOS+B,R200 t 

BD+28 4211 

20020904 

52521.83 

-2.4 

3200-10000 

CA+CAFOS+B,R200 t 

BD+28 4211 

20020905 

52522.82 

-1.4 

3200-10000 

CA+CAFOS+B,R200 t 

BD+28 4211 

20020906 

52523.85 

-0.3 

3500-10000 

CA+CAFOS+B,R200 

BD+28 4211 

20020907 

52524.89 

+0.7 

6200-10000 

CA+CAFOS+B,R200* 

BD+28 4211 

20020908 

52525.88 

+ 1.7 

3700-9300 

CA+CAFOS+B,R200 t 

BD+28 4211 

20020910 

52527.81 

+3.6 

3200-10000 

CA+CAFOS+B,R200 t 

BD+28 4211 

20020911 

52528.81 

+4.6 

3450-7600 

A1.82+AFOSC+gm4 


20020912 

52529.88 

+5.7 

3250-9800 

CA+CAFOS+B,R20 t 

BD+28 4211 

20020916 

52533.84 

+9.6 

3300-8700 

INT+IDS+R150V 

Kopff 27 

20020918 

52535.03 

+10.8 

3450-9600 

D1.54+DFOSC+gm4 

LTT 7987 

20020919 

52536.87 

+12.7 

3450-8800 

INT+IDS+R150V f 

Kopff 27 

20020922 

52540.00 

+15.8 

3300-9750 

D1.54+DFOSC+gm5 

LTT 7987 

20020923 

52540.86 

+16.7 

3300-9700 

INT+IDS+R150V 

Kopff 27 

20020924 

52542.00 

+17.8 

9400-24900 

NTT+SofI 

HIP 95550 

20020926 

52544.02 

+19.8 

3500-8800 

D1.54+DFOSC+gm5 

LTT 7987 

20021001 

52548.75 

+24.6 

3500-7500 

A1.82+AFOSC+gm4 

BD+28 4211 

20021010 

52557.84 

+33.6 

3600-9300 

INT+IDS+R150V 

Kopff 27 

20030410 

52739.17 

+215.0 

3500-8000 

TNG+LRB 

Feige 66 


Notes. * Epoch relative to B max , which occurred on MJD=52524.2 (Pignata et al. 2004). A1.82 = Asiago 1.82m telescope, CA = Calar Alto 
2.2 m telescope. INT = Isaac Newton Telescope. D1.54 = Danish 1.54 m telescope, NTT = ESO New Technology Telescope. TNG = Telescopio 
Nazionale Galileo. For the —8.3 d spectrum we used the flux standard taken with the same set-up on the previous night, while for the +4.6d 
spectrum we used a sensitivity function appropriate for the setup, kindly provided by N. Elias de la Rosa, f data not flatfielded; J missed target 
with blue grism. 


photospheric evolution, the differences become more apparent 
with time. 

3.2. Carbon in pre-maximum spectra 

Carbon is a diagnostically interesting element in pre-maximum 
spectra. Its abundance allows us to place limits on the compo¬ 
sition of the progenitor white dwarf. In addition, the amount 
of carbon constrains the extent of burning in the outer lay¬ 
ers thereby testing various explosion scenarios. For instance, 
current delayed-detonation models have very little carbon be¬ 
low about 30000 km s -1 , while the 1-D deflagration model 
(W7) has unburnt carbon present down to velocities of about 
14000 kms _1 (e.g. Branch et al. 2003). 

The first reported detection of carbon in a SN la was in the 
—14 d spectrum of SN 1990N (Fisher et al. 1997). However, 
carbon lines have also been identified at —10 d in the spectra of 
SN 1994D and at -9d in SN 1998aq (Branch et al. 2003). The 
strongest line is at A6580 A with weaker lines at A4267 A and 
A7234 A. In Fig. 4, we show tentative evidence for the pres¬ 
ence of the CII A6580 A line present up to about a week before 
maximum light at about 13000 kms -1 . There is no evidence 


of the CII A7234 A line (Fig. 4b) but this would, in any case, 
be expected to be much weaker. 

3.3. Expansion velocities 

As the photosphere recedes into deeper and slower moving ma¬ 
terial, the minima of the absorption features shift redwards with 
time. As already mentioned, this can provide a powerful “tomo¬ 
graphic” picture of the SN ejecta. However, very strong lines 
may be less useful in this regard. Patat et al. (1996) point out 
that the velocities inferred from the strong features (e.g. the 
characteristic Sill line) are only representative of the photo- 
spheric velocity at early phases when the amount of material 
above the photosphere is small; as the photosphere retreats, it 
encounters flatter density gradients which results in strong lines 
being formed over a large velocity range, thus overestimating 
the photospheric velocity. In contrast, weak lines can provide 
good indicators of the photospheric position. 

In SN 2002er, we examined the velocity evolution of the 
strong Sill A6355A and Call H&K lines and the weak SII 
A5640A line. The results are shown in Fig. 5 together with 
those of other supernovae for comparison. All velocities have 
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X(A) X(A) 

Fig. 4. Expanded view of several epochs of pre-maximum spectra 
around the CII A6580.7234A lines. The 02er spectra have been 
smoothed and shifted vertically for clarity, (a) A short vertical dash 
shows the position of the weak feature seen in the earliest spectra, 
(b) The dash shows the location of the A7234A line at the veloc¬ 
ity (13000 kms" 1 ) inferred from panel (a). The earliest spectra of 
SN 1994D and SN 1998aq are shown for comparison. 

been corrected for the recession velocities of the respective 
host galaxies. As expected, the weak SII line exhibits lower 
velocities than do stronger lines. The strong Sill feature in 
the earliest spectrum of SN 2002er (—11.3 d), has its mini¬ 
mum centred at about 6120 A, with its blue edge lying at about 
22000 kms -1 . As has already been pointed out by several 
authors (e.g. Mazzali et al. 1993) such velocities are well in 
excess of the 15000 kms -1 limit in the unmixed W7 model 
(Nomoto et al. 1984). 

Lentz et al. (2000) have computed the effects of varying 
the metallicity in the C+O layer of the standard deflagration 
model, W7, on the emergent pre-maximum spectra of type la 
supernovae. They find that the strength, profile, and velocity of 
the Sill A6355 A feature is a function of metallicity. Their re¬ 
sults are plotted in Fig. 5. Clearly the range of velocities seen 
between a number of SNe la may indeed be due to metallic¬ 
ity variations. However, even with implausibly large variations 
in the C+O layer metallicity (xl/10 to xlO solar) we can¬ 
not reproduce the high velocities exhibited by SN 2002bo and 
SN 1984A. These large velocities could be due to other factors 
such as higher explosion energy. 

4. Synthetic spectra 

Comparison with observations of synthetic spectra based on 
specific explosion models provides a powerful tool for probing 
the physics of supernovae. Here we compare the SN 2002er 
early-time spectra with model spectra generated by a Monte 
Carlo (MC) code. The code is based on the Sobolev approx¬ 
imation, and it was discussed in a number of papers that the 
reader can turn to for details (Abbott & Lucy 1985; Mazzali & 



-10 0 10 20 
Epoch (d) 


Fig. 5. The velocity evolution of the Sill A6355 A, Ca H & K, and 
SII A5640 A features compared with other SNe as measured from the 
minima of the respective absorption features. The solid curves shown 
in the top panel are models of different metallicity (top: x 10 middle: 
xl bottom:xl/10 solar) obtained from Lentz et al. (2000). Typical 
error bars for SN 2002er are shown in the bottom left-hand corner of 
each panel. 


Lucy 1993; Lucy 1999; Mazzali 2000). The input parameters 
are the luminosity L, the epoch t exp where explosion occurs at 
texp = 0d , the photospheric velocity v p h, and a density distri¬ 
bution from an explosion model. The code iteratively computes 
the radiation field and the level populations in the SN ejecta 
under the assumption that there is no net exchange of energy 
between matter and radiation above a sharply defined photo¬ 
sphere. This means that no deposition of fast particles created 
in radioactive decays occurs above the photosphere, and that 
absorptive continuum opacities are ignored. This is a reason¬ 
able approximation at the earliest phases of a SN la. 
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The luminosity is chosen to match the observed value and 
the velocity adjusted to obtain a good match to the observed 
line blueshifts. The epoch relative to explosion is determined, 
inevitably with some uncertainty, from the time of observed 
maximum. This is because first observations are usually at least 
a few days after the SN exploded. Different models can be 
chosen for the density distribution. In this paper we base our 
calculations on W7(Nomoto et al. 1985), which is known to 
give a good description of typical type la SNe. Energy pack¬ 
ets are released at the photosphere and travel through the SN 
envelope, where they can undergo line absorption and reemis¬ 
sion (branching) and electron scattering. When a converged 
temperature structure is established and appropriate excitation 
and ionisation conditions are computed (using a modified neb¬ 
ular approximation (Abbott & Lucy 1985; Mazzali & Lucy 
1993)), a synthetic spectrum is obtained using a formal inte¬ 
gral scheme. 

We calculated two spectral models for SN 2002er, one for 
an epoch near maximum light (day -2.4) and one for the earli¬ 
est epoch at-11.3 days. These two epochs were chosen to cover 
a period in which the physical properties change significantly. 

In the first iteration of the day -2.4 spectral model we 
adopted a typical value for the photospheric velocity at this 
epoch, of 9500 kms -1 . The distance modulus (/i = 32.9) and 
the reddening ( E{B — V) = 0.36) were taken from Pignata et 
al. (2004). This initial set of parameters provided a poor match 
to the observed spectrum. The problem was that, in order to 
match the observed bolometric luminosity of logioT ~ 43.2 
[erg s -1 ] (Pignata et al. 2004) and the observed line velocity 
(this work) the model temperature near the photosphere had the 
rather high value of Tbb ~ 14, 000 K. Consequently the ioni¬ 
sation was much too high. For example, strong Si III lines were 
predicted but not observed. The large observed line ratio of Si II 
5970 A to Si II 6350 A also indicates a rather low temperature. 
Nugent et al. (1995) derived a temperature of ~ 9000 K for 
SNe la with a line ratio similar to that of SN 2002er. 

In order to try to resolve this discrepancy between the 
model and observations, we explored the range of parameter 
space available through the uncertainties in reddening, distance 
and explosion epoch. We found that even if we set the red¬ 
dening at the 1 a lower limit of E{B — V) = 0.31 (Pignata 
et al. 2004), we make little difference to the discrepancy. The 
best match was obtained by pushing the distance modulus to 
its lcr lower limit of /j = 32.7, while keeping the formal 
value of the reddening, E(B — V) = 0.36, and using a rise¬ 
time of 20.Od. This risetime is 1.3d longer than estimated by 
Pignata et al. (2004), but is closer to the typical SN la rise¬ 
time of f,B max = 19.5 d (Riess et al. 1999). This departure from 
the Pignata et al. (2004) value can be justified by the fact that 
the photometric data are not available earlier than -7.3 days. 
Nevertheless, even with these more extreme parameter values, 
the temperature was still found to be rather too high. The best 
match parameters are login/, = 43.15[ergs -1 ] and a velocity 
at photophere of 9200kms -1 , leading to a photospheric tem¬ 
perature of Tbb = 12,460 K. The model thus obtained is com¬ 
pared with the observed spectrum in Fig. 6b. 

The main features are labelled with their rest wavelength 
in Fig. 6. We identify strong Si II lines at 6360 A, 5970 A, and 




Fig. 6. Observed spectrum and corresponding model of SN 2002er at 
using E(B — V) = 0.36 and fj, = 32.70 with / exp = 8.7 d (top panel) 
and /exp = 17.6 d (bottom panel). 


at 4130 A. Additionally, the deep absorption near 4800 A con¬ 
tains a significant amount of Si II. The rest of this feature is 
caused by Fe II. Weak Si III lines can be seen in the observed 
spectrum at 4570 A and 5740 A, but they are both very strong 
in the model due to the high temperature. S II produces the dou¬ 
ble line feature near 5400 A. The absorption at ~ 4400 A is due 
to Mg II 4481 A, as well as the one near 7500 A. Finally we 
recognise the strong Ca II H&K lines and the Ca II IR triplet. 
Although it is not as pronounced as in the case of SN 2002bo 
(Stehle et al. 2004), a high-velocity component in the Ca II lines 
can be seen that is not reproduced by the model, which uses the 
density distribution of W7. 

In constructing the spectral model for -11.3 days, we could 
not obtain the flux directly from observation since the earliest 
photometry epoch was at -7.3 days. [One possibility is to ex¬ 
trapolate the light curve to the early epochs using a quadratic 
function in time (Riess et al. 1999) (L oc t^ xp ). This procedure 
gives logioT — 42.75[ergs -1 ]. This is however too bright, 
leading to a very hot model. An alternative approach is to used 
the synthetic light curve fit of Pignata et al. (2004) to extrapo¬ 
late backto-11.3 days. This yielded logioT ~ 42.6[ergs -1 ]. 
This allowed the model to produce much more reasonable re¬ 
sults. The best match was obtained with a velocity at the pho¬ 
tosphere of 13000kms -1 , leading to Tbb = 10950K. The 
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model spectrum is compared with the observed spectrum in 
Fig. 6. Although the agreement is good, it should be kept in 
mind that the synthetic light curve shape may change when a 
more realistic element distribution is adopted compared with 
that used in Pignata et al. (2004). It is interesting that in 
SN 2002er, as in SN 2002bo, the temperature rises as we ap¬ 
proach maximum light. This behaviour is not typical of all 
SNe la (Benetti et al. 2004b) 

While the strong Si II line absorptions are well reproduced. 
Sill 5970A is too weak owing to the temperature still being 
too high. The S II double line feature near 5400 A is not well 
reproduced in its blue part. Since these lines come from the 
same lower levels, there may be a problem in the model line 
list, in particular the gf values. The line at 4700 A is not 
reproduced by the model. This feature was detected in other 
SNe la and can possibly be assigned to a high-velocity absorp¬ 
tion of Fe II (Hatano et al. 1999). All the other lines are well 
reproduced and have already been discussed for the day -2.4 
model. 

The above analysis indicates shortcomings in the mod¬ 
elling procedure. Specifically, owing to the homogenous abun¬ 
dance distribution of the model, to reproduce the strong Si II 
absorptions we had to introduce a rather high Si abundance 
(63% by mass) throughout the envelope. But due to the high 
photospheric temperature, the Si is almost completely doubly- 
ionised near the photosphere at maximum, severely overpro¬ 
ducing the Si III features. Only further out in the envelope, 
where the temperature is sufficiently low, are the required 
Si II absorptions formed. A similar problem was encountered 
with SN 2002bo (Benetti et al. 2004a). However, Stehle et al. 
(2004) find that, in a more recent analysis of SN 2002bo by 
the “Abundance Tomography” method, it is possible to reduce 
the Si abundance near the photosphere and yet simultaneously 
keep a desirably high Si abundance at the outer radii Stehle et 
al. (2004). Furthermore, Stehle et al. (2004) also find that im¬ 
proving the description of the radiation field at the inner bound¬ 
ary reduces the overestimate of the emergent flux in the red, 
leading to a decrease of the total luminosity and hence a lower 
temperature near the photosphere. In the case of SN 2002er, 
however, we also note that if the reddening and/or the distance 
are reduced by more than the 1 a error boundaries estimated by 
(Pignata et al. 2004) then greater consistency with the spectral 
models is achieved. At present, there is no independent evi¬ 
dence that such reduction should be applied. 

We conclude that the maximum light spectrum of 
SN 2002er can be fairly well modelled using a homogeneous 
abundance distribution in the atmosphere provided we push the 
adopted distance and risetime close to the observational lim¬ 
its. Modelling of the earliest spectrum depends sensitively on 
the adopted flux level, but unfortunately observations do not 
directly provide this value. However, extrapolation by means 
of a bolometric light curve model fit appears to provide a rea¬ 
sonable estimate of this quantity. Future progress will require 
a more realistic model which includes a stratified abundance 
distribution. 



Fig. 7. Near-infrared spectrum of SN 2002er at +18 d compared 
with spectra of other SNe la at comparable epochs. The spectra of 
SNe 1998bu, 1996X, 1999ee, and 2002bo have been taken from Jha et 
al. (1999). Hernandez et al. (2000), Hamuy et al. (2002), and Benetti et 
al. (2004a) respectively. The spectra are shown in the respective rest- 
frames of the supernovae and have been de-reddened using the values 
in Table 3. The epochs relative to maximum light are shown next to the 
name and the Amis values are noted just below. Prominent features 
are indicated following Marion et al. (2003). 

5. Near-IR spectroscopy of SN 2002er 

Early-time IR spectroscopic observations can provide a par¬ 
ticularly powerful diagnostic of the dynamic model because 
they probe different depths at the same epoch within the 
exploded white dwarf through the strongly variable line¬ 
blanketing opacity (Wheeler et al. 1998). Furthermore, the 
large number of lines in the optical/UV regions, coupled with 
the strong Doppler-broadening (several thousand kms -1 ) can 
make species identification difficult. Strong extinction effects 
often introduce additional error. In contrast, the near-IR has 
fewer lines and reduced sensitivity to extinction uncertainty, 
allowing firm line identification and accurate measurement of 
line strength and evolution. In addition, since there are much 
fewer strong resonance lines, it means that, in general, for a 
line to be detectable it needs a significant element abundance. 
In other words, IR lines suffer less from contamination by pri¬ 
mordial elements. 

After maximum light there is a rich mixture of Fe-group 
lines in the IR spectrum, allowing us to explore the innermost 
regions of the supernova and to ascertain the extent of incom¬ 
plete silicon burning. Wheeler et al. (1998) find that the post- 
maximum-light features at ~ 1.5—1.7 and 2.2-2.6 //m are pri¬ 
marily due to iron group elements of the second ionisation 
stage. The blue extent of the feature at ~ 1.55 /mi defines the 
transition from partial to complete Si burning (Marion et al. 
2003). Due to the high optical depth, this feature is formed 
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Fig. 8. The nebular spectrum of SN 2002er compared with the few 
spectra that are available of other SNe la at comparable epochs. The 
epochs relative to maximum light are indicated next to the names 
and the Amis values are indicated just below. The spectra have been 
corrected for the redshift of their respective host galaxies and have 
been dereddened in the same manner as for Fig. 3. F\ is in units of 
KT 16 erg cm -1 s -1 A -1 . The SN 1998bu, SN 1990N, and SN 1996X 
spectra have been shifted by +5.2, +4.0, and +3.2 respectively. The 
SN 1996X spectrum has been binned by a factor of 10 in an attempt 
to compensate for the very low signal in the 6000-7000 A region. 


near the outer boundary of the Fe/Co/Ni core and demarcates 
the the Fe/Si edge (Fig. 7) Wheeler et al. (1998) particularly 
note that the minimum at ~1.7 /xm provides a diagnostic tool 
for assessing the velocity spread of the region containing ra¬ 
dioactive materials which, in turn, depends on the underlying 
explosion scenario. The 2.2-2.6 /xm features also contain sig¬ 
nificant contributions from intermediate mass elements such as 
Si. Indeed, these K-band features are formed in the transition 
region between complete and incomplete silicon burning and 
are sensitive to the composition structure and hence the explo¬ 
sion model. 

Due to a combination of poor conditions and instrument 
availability, only a single epoch of near-IR spectroscopy was 
obtained, on 2002 September 24 (+17.8 d) using the NTT/SofI 
low-resolution blue and red grisms. The spectrum was re¬ 
duced using standard procedures in the FIGARO 4 environ¬ 
ment. Wavelength calibration of the 0.95-1.64 /xm spectrum 
was carried out using a Xe arc; as no arc frame was avail¬ 
able for the 1.53-2.52 /xm spectrum, we used the OH sky lines. 
Flux calibration was carried out with respect to the bright F6V 
star HIP 95550 and compared with near-IR photometry taken 
at +19.9 d (Pignata et al. 2004). The spectrum, scaled by 20% 
to match the photometry, is shown in Fig. 7. For comparison, 
we also show the rest-frame near-IR spectra of other SNe la at 
comparable epochs. Owing to the sparsity of near-IR spectra 


at the appropriate epoch (~3 weeks), we include two spectra 
of SN 2002bo that bracket the phase of the SN 2002er near-IR 
observation. 

Inspection of Fig. 7 shows that the IR spectrum of 
SN 2002er is fairly typical for the epoch considered. As an¬ 
ticipated for this epoch, the spectrum is dominated by the in¬ 
creasingly exposed Fe-group elements. The 1,7/xm feature is 
considerably shallower in SN 2002er than in SN 1996X and 
is more reminiscent of the morphology seen in the somewhat 
faster decliners such as SN 1998bu and perhaps, SN 1999ee. 
As indicated above, this may be interpreted as a larger-than- 
average velocity spread in the radioactive materials. However, 
the measured Doppler blueshift in the 1.5/xm edge (rest wave¬ 
length 1.57^tm, Marion et al. 2003) is 8300±500 kms -1 for 
SN 2002er compared with 10500±500 kms" 1 for SN 1996X 
at +19d. This suggests that silicon burning actually extended to 
higher velocities in SN 1996X. 

6. The nebular spectrum at +215 days 

In a Type la supernova, by about 60-100 days post-explosion 
the photosphere has receded to the centre, and the ejecta 
have become transparent to the optical/IR emission i.e. the 
forbidden-line-dominated nebular phase has been reached. The 
transparency of the ejecta means that the physical condi¬ 
tions and constituents of the innermost regions can be directly 
probed. Furthermore homologous expansion means that line 
profiles provide valuable information about the distribution of 
the nucleosynthesised materials. In particular, the mass and ve¬ 
locity distribution of 56 Ni and its products can provide infor¬ 
mation on the nature of the explosion. 

We obtained a single nebular phase spectrum at +215 d post 
/7-maximum. The spectrum was reduced as described in Sec. 2 
and is shown in Fig. 8 together with spectra of other type la su¬ 
pernova at comparable epochs. As with the earlier epoch spec¬ 
tra, the form of the late time spectrum of SN 2002er is not un¬ 
usual. The spectrum is dominated by forbidden lines of singly 
and doubly ionised Co and Fe. The FWHM of the 4700 A fea¬ 
ture (~13000 km s -1 ) is consistent with the relation in Mazzali 
et al. (1998) that links the width of nebular lines to the lumi¬ 
nosity of the supernova. 

In order to exploit the information available in the nebu¬ 
lar spectrum, we compared it with the non-LTE nebular spec¬ 
tral model described in Bowers et al. (1997), to which we re¬ 
fer the reader for details. Briefly, the model comprises a uni¬ 
form density, homologous expanding sphere containing iron, 
cobalt, and sulphur with the relative abundances as specified 
by radioactive decay. Forbidden lines of these elements emit¬ 
ted by singly and doubly-ionised species are predicted to dom¬ 
inate the nebular phase at optical and near-IR wavelengths. 
Consequently only these two ionisation stages were included 
in the model. The distance and reddening were set at /i = 32.9 
and E{B — V ) = 0.36 respectively (see above). The 56 Ni mass, 
expansion velocity, ionisation degree and temperature are free 
parameters. The best-match model is shown in Fig. 9, com¬ 
pared with the observed spectrum. The model parameter values 
are T = 6600 K, expansion velocity ~90001< m s “ 1 fractional 
abundance of SII=0.5, and M( 56 Ni)= 0.69 M©. These clearly 
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yield a fair match to the observed spectrum. We also note that 
the derived mass of 56 Ni is close to the value obtained from the 
bolometric light curve (0.7 M 0 , Pignata et al. 2004). We con¬ 
clude that (i) the masses of 56 Ni derived by the two methods 
are consistent with each other, and (ii) the fraction of neutral 
and triply-ionised iron-group elements is indeed quite small. 

7. Summary 

We have presented extensive optical spectroscopy of 
SN 2002er at epochs spanning —11 to +215 days. The 
spectroscopic coverage at epochs up to +34 days is especially 
noteworthy. In most respects, SN 2002er bears the hallmarks 
of being a typical, albeit significantly reddened, type la super¬ 
nova. Its spectral evolution was similar to those of other normal 
type la SNe such as SN 1994D and SN 1996X that suffered 
from minimal reddening. Nevertheless, differences are seen 
between the coeval spectra of SN 2002er and other normal 
SNe la. These are most pronounced in the pre-maximum 
spectra. 

We compared the spectra at two early-time epochs with a 
homogeneous spectral synthesis model based on the deflagra¬ 
tion model W7. The best match was obtained only by push¬ 
ing the distance modulus to its la lower limit of // = 32.7, 
and using a risetime as high as 20.Od. Even then the model 
temperature is rather too high. A similar problem was encoun¬ 
tered at both epochs. Nevertheless, most of the features of the 
early-time spectra were quite successfully reproduced by the 
model. Future modelling improvements will include stratified 
abundances, plus a better description of the radiation field. 

Only one early-time IR spectrum was obtained. As with the 
optical spectra, this showed that SN 2002er was fairly typical 
but with some differences from other events. 

An optical spectrum obtained at +215 d spectrum was mod¬ 
elled using a simple one-zone nebular code. This provided a 
fair match to the observed spectrum. A o6 Ni mass of 0.69 M q 
was inferred - consistent with that derived from the light curve 
(O.7M0 Pignata et al. 2004), and indicating that the fraction of 
neutral and triply-ionised iron-group elements is quite small. 

The principal contribution of this paper to the study of 
the SN la phenomenon has been its provision of an exten¬ 
sive, frequently-sampled set of early-time spectra. These spec¬ 
tra span the crucial phase during which the photosphere re¬ 
cedes from the outermost unburned layers into the innermost 
iron-group regions of the core. This provides a superb dataset 
against which explosion scenarios will be tested via “tomo¬ 
graphic” spectral synthesis models. Such studies will be vital 
in order to understand the significance of the small differences 
seen between different type la supernovae. 
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Fig. 2. Spectral evolution of SN 2002er during the photospheric phase. The spectra shown here have neither been corrected for the redshift 
of UGC 10743, nor have they been corrected for reddening. The spectra have been shifted vertically for clarity; these shifts are indicated in 
parenthesis next to the epoch of the spectra (see also Table 2). Note that the quality of the spectra beyond ~9000 A is poor. The prominent 
telluric features are marked by earth symbols. 
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Fig. 3. Spectra of SN 2002er compared with other type la SNe at four representative epochs. The spectra have been shifted vertically by arbitrary 
amounts and have been corrected for the respective redshifts of the host galaxies. All spectra (except those of SN 1998aq) have been corrected 
for extinction using the values in Table 3 and the relation of Cardelli et al. (1989). The Amis values shown in panel (b) have been corrected 







